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ABSTRACT 

The MAGIC telescope took data of very high energy 7-ray emission from the blazar Markarian 421 
(Mkn 421) between November 2004 and April 2005. We present a combined analysis of data samples 
recorded under different observational conditions, down to 7-ray energies of 100 GeV. The flux was 
found to vary between 0.5 - 2 Crab units (integrated above 200 GeV), considered a low state when 
compared to known data. Although the flux varied on a day-by-day basis, no short-term variability was 
observed, although there is some indication that not all nights are in an equally quiescent state. The 
results at higher energies were found to be consistent with previous observations. A clear correlation 
is observed between 7-rays and X-rays fluxes, whereas no significant correlation between 7-rays and 
optical data is seen. The spectral energy distribution between 100 GeV and 3 TeV shows a clear 
deviation from a power law, more clearly and at lower flux than previous observations at higher 
energies. The deviation persists after correcting for the effect of attenuation by the extragalactic 
background light, and most likely is source-inherent. There is a rather clear indication of an inverse 
Compton peak around 100 GeV. The spectral energy distribution of Mkn 421 can be fitted by a one- 
zone synchrotron self-compton model suggesting once again a leptonic origin of the very high energy 
7-ray emission from this blazar. 

Subject headings: gamma rays: observations, BL Lacertae objects: individual (Mkn 421) 
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1. INTRODUCTION 

Mkn 421 (redshift z = 0.030) is the closest known 
and, along with Mkn 501, the best studied TeV 7-ray 
emitting blazar. It was the first extragalactic source 
detected in the TeV energy range using imaging at- 
mospheric Cherenkov te lescopes (IACTs) ([Punch et al.1 
Il992t iPetrv et al.l I1996D . Mkn 421 is currently the 
source with the fastest observed flux variations among 
TeV 7-ray emitters. So far it has shown flux varia- 
tions larger than one order of magnitude, a nd occasional 
flux doubling times as shor t as 15 min (|Gaidos et al.l 
Il996t lAharonian et al.l l2002h . Variations in the hard- 
ness of the TeV 7-ray spectru m during flares were re - 
ported by severa l group s (e.g. iKrennrich et alj ([2002) ; 
lAharonian et al.l ((2005)). Simultaneous observations 
in the X-ray and GeV- TeV bands show strong ev- 
idence for flux correl ation (jKrawczvnski et ail 120011 : 
IBlazeiowski et al.ll2005l ). 

u Universitat de Barcelona, E-08028 Barcelona, Spain 
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Mkn 421 has been detected and studied in all acces- 
sible wavelengths of the electromagnetic spectrum from 
radio waves to very high energy (VHE) 7-rays. The over- 
all spectral energy distribution (SED) shows a typical 
two bump structure with the first peak in the keV en- 
ergy range and the second maximum at GeV-TeV ener- 
gies. The SED is commonly interpreted as beamed, non- 
thermal emission of synchrotron and inverse-Compton 
radiation from ultrarelativistic electrons, accelerated 
by shocks moving along the jets at relativistic bulk 
speed. Simple one-zon e synch r otron-self-Compton (SSC) 
model s (e.g. ICoppil (jl992t ); ICostamante fc Ghisellinil 
(2002)) describe the observation al results satisfactorily 
well. However, had ronic models (jMannheim et al]|1996t 
iMucke et al.l [2003) can explain the observed features, 
too. A way to distinguish between the different emission 
models is to determine the position of the second peak in 
the SED, using simultaneous time-resolved data over a 
broad energy range through multiwavelegth campaigns. 
This requires providing data in the as yet unexplored gap 
in the SED. 

The MAGIC telescope (Major A tmosph e ric G amma 
Imaging Cherenkov telescope; see iLorend ((2004)) , lo- 
cated on the Canary Island La Palma (2200 m asl, 
28°45'N, 17°54'W), completed its commissioning phase 
in early fall 2004. MAGIC is currently the largest IACT, 
with a 17 m diameter tessellated reflector dish consisting 
of 964 0.5 x 0.5 m 2 diamond-milled aluminium mirrors. 
Together with the current configuration of the MAGIC 
camera with the tri gger region of 2.0 degrees diameter 
(|Cortina et al.|[2005f ). this results in a trigger collection 
area for 7-rays of the order of 10 5 m 2 , increasing with 
the zenith angle of the observation. Presently the acces- 
sible trigger energy range spans from 50-60 GeV (at small 
zenith angles) up to tens of TeV. The MAGIC telescope 
is focused to 10 km distance - the most likely position 
for a 50 GeV air shower. The accuracy in reconstructing 
the direction of incoming 7-rays on an event by event 
basis, hereafter 7 point spread function (PSF), is about 
0.1 degrees, slightly depending on the analysis. 

The first physics observations in winter 2004/05 and 
in spring 2005 included observations of the well estab- 
lished TeV blazar Mkn 421. In total, 19 nights of data 
were taken on this source, the observation times per night 
ranging from 30 minutes up to 4 hours. Here we present 
the results from these observations, covering the energy 
range from 100 GeV to several TeV. We first describe 
the data set and analysis techniques in section O In sec- 
tion 02 we present the results and, finally, in section [4] 
we compare our results with other observations and in- 
terpret them in terms of different models. 

2. OBSERVATIONS AND DATA ANALYSIS 

The Mkn 421 data were taken between November 2004 
and April 2005, and divided into four samples, for rea- 
sons given below. Data taken before and after February 
2005 were treated separately, due to changes in the tele- 
scope hardware. Most of the data were taken at small 
zenith angles (ZA < 30°), i.e. at a low trigger energy 
threshold. However, observations made during 1.5 hours 
in a common campaign with the H.E.S.S. telescope sys- 
tem azm et al.l 120050 in December 2004 were taken 
at 42° < ZA < 55°. There were also different obser- 
vational modes: the standard mode for MAGIC is the 



ON-OFF mode, with equal time given to tracking the 
source in the center of the camera (ON), and tracking a 
sky region near the source but with the source outside the 
field of view (OFF). This provides a robust estimate of 
the background. In our observations, we considered the 
7-ray signal from Mkn 421 to be strong enough to obvi- 
ate OFF observations, and we estimated the background 
level from the ON data (see below). In Apr il 2005, part 
of the data were taken in the wobble mode (jDaum et al.l 
Il997f h In this mode, two sky directions, opposite and 
0.4° off the source, were tracked alternately for 20 min- 
utes each, which provides a simultaneous measurement 
of signal and background. In the wobble mode there is 
a priori no need for additional OFF data. 
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Fig. 1. — Alpha distribution for the combined data samples 
I+II+iII with E t i lrea f l = 150 GeV. A vertical line indicates the 
Alpha cut used to extract excess events. The black parabola 
is a fit to the Alpha distribution between 30 and 90 degrees 
and is used to estimate the background level between and 10 
degrees. 



The observation criteria and some important parame- 
ters of the four data samples are summarized in Table [TJ 
For each data sample a separate Monte-Carlo (MC) set 
of 7 events was sim u lated (CORSIKA versio n 6.023, 
iKnapp k. Heckl (|2004h ; iMaiumdar et all (|2005h ). taking 
into account the zenith angle of observation, the obser- 
vational mode, and the hardware setup of the telescope. 
The full data set corresponds to 29.0 hours ON-source 
observation time. Runs with problems in the hardware 
or with unusual trigger rates were rejected in order to 
ensure a stable performance and good atmospheric con- 
ditions. After removing these runs, the remaining obser- 
vation time was 25.6 h. 

For calibration, image cleaning, cut optimization, and 
energy reconstruction t he standard a nalysis techniques of 
the MAGIC telescope (iBreti [20051: iWaener et all 120051: 
lOaug et all [20051 were applied as shortly described be- 
low. The calibration of the raw data from the MAGIC 
camera uses a system consisting of fast and powerful 
LED pulsers emitting at three different wavelengths with 
variable light intensity. Absolute calibration is obtained 
by comparing the signal of the pixels with the one ob- 
tained from a carefully calibrated PIN diode, and is cross- 
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TABLE 1 

Results of the Mkn 421 data using the Alpha approach (see text for details). Samples I+II were recorded in November 
2004 - January 2005, while samples III+IV were taken in April 2005. 



sample 


on time 


ZA range [°] 


mode 


B thr [GeV] 


N on 


N off 


N cxceas 


sigma 


i 


4.63 h 


9.3 - 31.2 


ON 


150 


3761 


1878 ± 32 


1883 ± 69 


29.3 


ii 


1.53 h 


42.4 - 55.0 


ON 


260 


1086 


674 ± 25 


413 ± 41 


10.1 


in 


9.30 h 


9.2 - 27.5 


ON 


150 


8083 


4360 ± 49 


3723 ± 102 


38.9 


IV 


10.12 h 


9.4 - 32.4 


wobble 


150 


7740 


4532 ± 67 


3208 ± 111 


29.1 
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Fig. 2. — 8 2 distribution for the combined data sam- 
ples III+IV with i?th res h = 150 GeV after background 
subtraction. The black line is the MC-7 expectation for 
a point source. 
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3. — Sky map of excess events in the re- 
Mkn 421 for samples III+IV using the Disp 
method ( Domingo-Santaman'a ct al. 2005). The black cross 
indicates the reconstructed source position. Note that the ver- 
tical scale is in units of [events / (0.05x0.05 deg 2 )]. 



checked by analysing muon rings. The time resolution 
of the read-out system has been measured to be about 
700 ps for Cherenkov light flashes of 10 photo-electrons 
(ph. el.) per pixel, reaching 200 ps at 100 ph. el. Calibra- 
tion events are taken at 50 Hz, interlaced with normal 



data, using an external calibration trigger. 

The calibrated images are cleaned using so-called tail 
cuts: pixels are retained only if their reconstructed 
charge signals are larger than 10 ph. el. ('core pixels') or 
if their charges are larger than 5 ph. el. and they have 
at least one neighboring core pixel. The c amera image s 
are then reduced to image parameters as in (Hillas 1985), 
adding parameters describing the intensity concentration 
and asymmetry. 




10" 

Energy [GeV] 

Fig. 4. — Differential energy distribution of the Crab Nebula 
data sample from 2004 ( black circles) and 20 05 (grey circles) 
as measured by MAGIC (Wagner ct al. 2005). A fit by power 
law with a changing photon index to the MAGIC 2004 data 
is shown by the dashed black line. The analytical form of the 
fit and the fit p arameters ar e listed in the inlay. Wipple data 
(solid bl ue line. IHillas et al.| [[ 1998)1 and HEGRA data (solid 
red line. [Aharonian ct al. (2004)) are shown for comparison. 



For 7/hadron separation a multidimensional classifica- 
tion technique based on the Rando m Forest (RF) method 
(jBreimanl 120011: iBock et al.l 120041 ) was used. The RF 
method uses training data (randomly chosen data events 
and Monte Carlo 7s, representing background and sig- 
nal) to find a set of classification trees in the space of 
image parameters. Multiple trees are combined to form 
a generalized predictor by taking the mean classification 
from all trees. The predictor, called hadronness, spans 
a range between and 1, and characterizes the event 
images being less or more hadron-like. 

In our analysis, classical image shape parameters like 
Width, Length, Dist and Size were used as input pa- 
rameters. The cuts in hadronness for the 7/hadron sepa- 
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ration were trained for each data set separately, and were 
then chosen such that the overall cut efficiency for MC 7 
events remained about 50%. The corresponding hadron 
suppression is about 90-99%, improving with increasing 
Size of the events. 

A critical variable not used in the RF classification tree 
is Alpha, the angle between the major image axis and the 
line connecting the center of gravity of the image with 
the source position in the camera plane. In stand-alone 
IACTs, Alpha is commonly used, after all previously 
noted cuts, to extract the 7 signal from the data, and 
to estimate the level of background. For a point source, 
the Alpha distribution of the 7-like events is expected 
to peak at low values of Alpha, whereas for background 
events the distribution should be flat or slowly varying 
with Alpha. 

In the case of our ON-mode data, the background re- 
maining after 7/hadron separation was estimated from 
the Alpha distribution by performing a second order 
polynomial fit (without linear term) in the range be- 
tween 30° and 90° where no contribution from 7 events 
is expected (see Fig. [1]). The signal was then deter- 
mined as the number of observed events in the range 
Alpha < Alphao exceeding the fit extrapolated to small 
Alpha, where Alphao is energy dependent and has a typ- 
ical value of 15°. The significan ce of an excess is then 
calculated according to Eq. 17 in lLi fc Mai (|1983f) . 

In the wobble mode, the ON (source) data are de- 
fined by calculating image parameters with respect to the 
source position, whereas OFF data are obtained from the 
same events but with image parameters calculated with 
respect to the position on the opposite side of the camera, 
the antisource position. In order to avoid an unwanted 
contribution from 7-events in the OFF sample and to 
guarantee the statistical independence between the ON 
and the OFF samples in the signal region, the following 
procedure is applied: events with Alpha source < Alpha® 
(with Alpha source calculated with respect to the source 
position) are excluded from the OFF sample, and events 
with Alpha an usource < Alpha® (with Alpha ant isource cal- 
culated with respect to the antisource position) are ex- 
cluded from the ON sample. This cut assures that the 
Alpha distributions for ON and OFF events are statis- 
tically independent for Alpha < Alpha®. The Alpha 
approach was used to determine the excess events for all 
four data sets (Table [1]) . 

As an alternative to this classical approach using 
Alpha, the so-called 9 2 approach can be applied, an 
approach more common for the analysis of data from 
a system of IACTs like HEGRA or H.E.S.S. The an- 
gle 9 denotes the angular distance between the ex- 
pected source position and the reconstructed origin of 
the initial 7-ray. Since for a single IACT the an- 
gle 9 can not be reconstructed d i rectly, the so - called 
Disp method (iFomin et al.l 11994 iLessard et all 120011 : 
iDomingo-Santamaria et al.ll2005D was used to determine 
the source position in the camera plane, using position- 
independent image shape parameters. The number of 
excess events is then determined as the difference be- 
tween the 9 2 distributions for the source and back- 
ground, respectively, similar to the Alpha approach. The 
background-subtracted 9 2 distribution for samples III 
and IV is shown in Fig. [2] The average background was 
estimated from the wobble data themselves, by exclud- 



ing the sector of the camera affected by the presence of 
the strong source. The solid line in Fig. indicates the 
expectation from MC-7 events for a point source. Com- 
puting 9 2 also permits to produce sky maps in which for 
every 7-ray candidate an origin in the sky is assigned (see 
Fig. [3]). Note that our signal analysis relies on the Alpha 
approach throughout. 

These conservative analysis methods are known to pro- 
duce reliable results for energies above 100 GeV. The en- 
ergy regime below 100 GeV will require additional stud- 
ies, in particular concerning the background rejection. 
Thus, for our analyzed sample the Size parameter (to- 
tal amount of light of the image and in first order pro- 
portional to the energy) was required to be above 150 
photoelectrons. 

The energy estimation was performed using again the 
Random Forest technique, based on the image parame- 
ters of a MC 7 sample. This sample is statistically inde- 
pendent of the one used for the training of the 7/hadron 
separation cuts. Prior to the training of the energy esti- 
mation, loose (high-acceptance) cuts in hadronness and 
Alpha were applied to avoid a possible bias caused by 
outlier 7 events. 

The energy thresholds of the individual analyses (as 
given in Table Q} are defined as the peak in the differ- 
ential energy distribution of the MC-7 events after all 
cuts. Our analyses showed that we were able to extract 
excess events with energies ~50 GeV lower than the cor- 
responding peak value. 

The spectrum of the number of excess events in bins 
of true energy is determined from the spectrum in the 
estimated energy by an unfolding procedure. This pro- 
cedure corrects for the finite energy resolution and for 
biases in the energy estimation. The unfolding program 
package used in MAGIC allows unfold ing with a vari- 
ety of methods (|Anvkevev et alj (|1991h ). which differ in 
the way regularization is implemented. Unfolding results 
are only accepted if the results from the different meth- 
ods are consistent with each other and if some criteria 
are satisfied concerning the regularization strength, the 
size of the noise component and the \ 2 value. The lat- 
ter is a measure of the agreement between the expected 
"measured" spectrum from the unfolded spectrum and 
the actually measured spectrum. The unfolding result 
presented in this analysi s was obtained with an iterative 
method, as described in iBerterol (|1989f) . 

To demostrate the quality of the applied analysis and 
the good agreement with previous measurements by 
other experiments, we show in Fig. [4] the differential en- 
ergy spectrum of the Crab Nebula data ("standard can- 
dle" of VHE 7-ray atronomy). The data were taken in 
2004 and 2005 with observation conditions and telescope 
perfomance similar to those of the Mkn 421 data. Ad- 
ditional publications describing details of the calibration 
methods and the data analysis are in preparation. 

3. RESULTS 

3.1. The signal 

During the entire observation period, Mkn 421 was 
found to be in a low flux state compared with existing 
data (around 1 Crab unit for a flux integrated above 
200 GeV), but resulting in a clear signal in all four data 
samples. Fig. [T] shows the Alpha distribution of the 7- 
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Fig. 5. — Night-by-night light curve for Mkn 421 from November 2004 to April 2005. Left panels: data from November 2004 to April 
2005. Right panels: expanded data for 6 nights in April 2005. Upper panel: MAGIC data, night average of Mkn 421 above 200 GeV using 
samples I+III+IV. Middle panel: corresponding day-by-day X-ray counts as observed by the RXTE/ASM. Lower panel: Light curve of the 
optical flux of Mkn 421 as measured by the KVA telescope. 



candidates of the combined samples I, II, and III with 
an energy threshold of ~150 GeV. An excess of about 
7000 events was found, which, for the given background, 
corresponds to an excess of more than 49 standard devia- 
tions. The number of excess events and the significances 
for the individual samples are summarized in Table [TJ 

Fig. [3] shows a sky map produced with the Disp 
method using samples III and IV. The reconstructed 
source position from the sky map (Fig. [3J indicated 
by the black cross) is centered at RA=+llh04'19", 
DEC=38°11'41". The systematic error of the telescope 
pointing is 2'. The 7 PSF is indicated by a white cir- 
cle in the left bottom corner. The observed extension 
of Mkn 421 is compatible with the MC expectation of a 
point source, which can also be seen in Fig. O 

3.2. The light curve 

The integral fluxes above 200 GeV, averaged over each 
night of observation, are shown in the upper panels of 
Fig.[U Significant variations of up to a factor of four over- 
all and up to a factor two in between successive nights 
can be seen. Since sample II has an energy threshold of 
260 GeV it is not shown on the light curve. The relatively 
high analysis energy threshold of 200 GeV applied for the 
light curve ensures that the results are independent of the 
actual trigger thresholds during each night. In the mid- 
dle panels of Fig. O the corresponding flux in the X-ray 
band as observed by the All-Sky-Monitor (ASM 23 ) on- 
board the RXTE satellite is shown. In the lower panels 
of Fig. [5] the optical data taken by the KVA telescope 24 
on La Palma are shown. Note that the contribution of 
the host galaxy (appr. 8.0 mJy) has been subtracted. 
While the X-ray data show a moderate variability within 
the observation period, the optical flux stays almost con- 
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stant. 

For the 6 nights in April (MJD 53465 to 53471), the 
light curve above 200 GeV is shown in Fig. [6] in bins 
of 10 minutes. We also added the background rates for 
each night in the same binning, in order to demonstrate 
that the small variations in the excess rates and the daily 
changes are not caused by detector effects or atmospheric 
transmission changes. The vertical lines indicate the time 
in each night at which the observation mode was changed 
from ON to Wobble. The mean integral flux per night R 
in units of 10 _9 cm~ 2 s _1 and the quality of the fit con- 
stant are shown in the panels. The horizontal dashed 
line corresponds to the integral flux of the Crab Nebula 
above 200 GeV. Combining the findings from the intra- 
night light curves we conclude that we did not find signif- 
icant short-term flux variability within individual nights, 
despite the high sensitivity of MAGIC for such a search. 
Some of the nights, however, are less compatible with 
a constant flux than others, which might be an indica- 
tion of some activity, albeit unstructured and difficult to 
quantify. On the other hand, we observe significant day- 
to-day variation by up to a factor of two, and differences 
up to a factor of four in the full sample. 

3.3. The energy spectrum 
3.3.1. The measured spectrum 

For the spectrum calculation, we combined the entire 
data set because the differences between the fluxes on 
individual nights are rather moderate (see Fig. [5]) . The 
resulting averaged differential energy spectrum is shown 
in Tableland in Fig. [7] by filled grey boxes. The energy 
spectra extend from around 100 GeV to several TeV. The 
last spectral point at 4.4 TeV is an 95% upper limit. The 
error bars shown are statistical only. Systematic errors 
are estimated to be 18% on the absolute energy scale, 
which correspond to 44% on the absolute flux level for a 
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Fig. 6. — Light curve for 6 nights in 2006 April in 10 minutes binning. Upper panels: flux above 200 GeV. Mean rate R in units 
of (10 _9 cm — 2 s _1 ) and the quality of the fit by a constant are shown in the panel. Lower panels: mean background rate R D g per 

minute after cuts. Note the rising background rate towards the end of each observation slot, related to the rising moon. R D g and the 
quality of the fit by a constant are shown. The vertical dotted lines indicate the time of the switchover from the ON observational 
mode to the Wobble mode. The dotted horizontal line indicates the Crab Nebula integral flux above 200 GeV as measured by 
MAGIC dWagner et al.ll2005h . 
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photon index of 2.2. The systematic error on the slope is 
estimated to be 0.2. The attenuation of the VHE photons 
by intergalactic low energy photons and the determina- 
tion of the intrinsic spectrum of Mkn 421 are discussed 
below. 

3.3.2. -y-ray absorption by the EBL 

The VHE photons from Mkn 421 cross -400 million 
light years on their way to Earth. They interact with the 
low energy photons o f the ex t ragalactic background ligh t 
(EBL. see (Nikishovl dl96l; iGould &, Schreden (fl966h : 
IStecker et all (|1992| ): lHauser fc Dwekl (|2001l )T consist- 
ing of redshifted star light of all epochs and reemis- 
sion of a part of this light by dust in galaxies. The 
most common reaction channel between VHE 7-rays and 
the low energy photons of the EBL is pair production 
7vhe + 7ebl — » e + e~, a reaction which has its largest 
cross section when the center of mass energy is roughly 
3.6 times larger than the threshold energy of 2m e c. The 
intrinsic (de-absorbed) photon spectrum, dN/dEi, of a 
blazar located at redshift z is given by: 

dN/dE l = dN/dE obs x exp[r 77 (S, z)], 

where dN/dE b s is the observed spectrum and r 77 (£', z) 
is the optical depth. The distance to Mkn 421 implies 
that t he optical depth (e.g. Eq. 2 in lDwek fc Krennrichl 
(2005)) strongly depends on the shape and absolute pho- 
ton density of the EBL between 1 and 30 /jm wavelength. 
A rather complicated distortion of the intrinsic spectrum 
takes place above ~100 GeV. Although the calculation of 
the optical depth is straightforward, the spectral energy 
distribution of the EBL is uncertain. Direct measure- 
ments of the EBL are difficult because of the strong fore- 
ground emission consisting of reflected sunlight and ther- 
mal emission from zodiacal dust par ticles. Hence, many 
measurements lead to upper limits (|Hauser et al.|[l99g : 
iDwek fc Arendtl Il998l ). Several measurements claimed 
a direct detec tion of the EBL, but some of them are 
contr oversial ijMatsumoto et alJ 120051 : iFinkbeiner et al.1 
2000) . An alternative method to determine the EBL are 
fluctuation analyses of the measured radiation. Since a 
part of the EBL originates from discrete sources, fluc- 
tuations in the number of sources in an observer's field 
of view will produce fluctuations inthe measured back- 
groun d (|Kashlinskv et al.l fl996: Kashlins kv fc Odenwaldl 
2000) . A third method is the galaxy number counting in 
the deep field surveys, which provides robust lower limits 
to the SEP of the EBL (lElbaz et alJl2002t iMetcalfe et all 
120031: iFazio et all l2004t iMadau fc Pozzettil I2000D . The 
results of these methods and measurements are summa- 
rized in Fig. [51 

In principle, upper limits on the EBL can also be deter- 
mined from observed 7-ray spectra from medium to high 
redshift TeV blazars. Under assumptions that the recon- 
structed TeV blazar spectrum is not too hard and it does 
not have a pile-up at high energies , the EBL level can 
be constrained (see lHauser fc Dwekl ([200l for summary 
and lAharonian et al.l (|2006l ) for latest results). However, 
since the measured spectrum of Mkn 421 is much softer 
that the one of Mkn 501, which is located at similar 
redshift, the softness seems to be intrinsic. In addition, 
the data in this paper extend up to 3 TeV only (his- 
torical data of Mkn 421 extend up to 20 TeV), which 
further weakens possible constraints from such a nearby 
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Fig. 7. — Differential energy distribution for Mkn 421 averaged 
over the whole data sample. The measured energy spectrum is 
shown by the grey full squares and the de-absorbed spectrum by 
the black full circles. The spectral point at the highest energy is a 
95% upper limit. The grey shaded area corresponds to a system- 
atic error from a slope error of ±0.2 as quoted in the text. The 
black solid line indicates the best fit to the de-absorbed spectrum 
by a power law with exponential cut-off; its parameters are listed 
in the inl ay. For comparison reasons, the measured Crab Nebula 
spectrum (Wagner et al,| [2005D is shown by the grey dotted line. 



TABLE 2 

Averaged measured differential energy spectrum of 
Mkn 421. derived from the data set presented in the 
Paper. The last point is an 95% upper limit. 



energy bin [GeV] 



differential flux dN/dE 



lower bin 


mean 


upper 


[photons/ (TeV cm 2 


»)] 


limit 


energy 


bin limit 


108 


134 


167 


(3.72 ±0.34) x 10~ 


-9 


167 


208 


259 


(1.21 ±0.04) x 10" 


-9 


259 


321 


402 


(3.77 ±0.15) x 10~ 


10 


402 


498 


623 


(1.32 ± 0.05) x 10- 


10 


623 


770 


965 


(3.63 ± 0.19) x 10- 


11 


965 


1192 


1497 


(8.95 ± 0.71) x 10- 


12 


1497 


1845 


2321 


(2.26 ± 0.27) x 10- 


12 


2321 


2856 


3598 


(2.88 ± 1.20) x 10- 


13 


3598 


4429 


5579 


< 1.10 x 10- 13 





source. We therefore do not try to constrain EBL using 
this Mkn 421 data set. 

Inst ead, we adopt the recent model of iPrimack et al.l 
(2005), scaled up by a factor 1.5 (which is within the 
model uncertainties), to match lower limits set by the 
Spitzer missi o n and ISOCAM in t he ran g e of 4 to 15 /xm 
(IFazio et al.l 120041; lElbaz et alJ l2002t iMetcalfe et ail 
l2003h . The resulting EBL spectrum is shown in Fig. [8] by 
the black curve. This EBL spectrum a g rees with alterna - 
tive models (e.g. iKneiske et al I (l200l; iPei et all (|1999T) : 
iBlain et alJ (|1999ft : IStecker et all (|2006t) 1 which are de- 
signed to predict the EBL today. It is also very close 
to the upper limits inferred fr om arguments on AGN 
spectra (|Aharonian et al.1 120061 ) . Using this EBL spec- 
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Fig. 8. — Energy density of the extragalactic background light (EBL). Direct measurements, galaxy counts, low and upper limits 
are shown by different symbols as described in the legend. The black solid curve is the EB L spectrum as in Primack et al. 120051) 
for z=0 but upsc aled by a factor 1.5 to match low limits derived from the galaxy counts (p Efbaz e t al. 2002; Metcalfe et alj 12*00*3 : 
Fazio et al. 2001). 



TABLE 3 

Systematic study of the fit parameters on the 
de- absorbed spectrum of mkn 421. the fitted fuction is a 
power law with an exponential cut-off: 
dN/dE = /V u (E/0.2TeV)- a exp(-E/E cutoS ). We show fit 

VALUES ON THE PHOTON INDEX, a, AND THE CUT-OFF ENERGY, 
-Ecutoffi FOR FOLLOWING ASSUMPTIONS: NOMINAL VALUES (A), A 
SYSTEMATIC SHIFT BY +18% IN THE VHE ENERGY SCALE (B), A 
SYSTEMATIC SHIFT BY -18% IN THE VHE ENERGY SCALE (C), A 
SYSTEMATIC SHIFT BY +18% IN THE VHE ENERGY SCALE AND, IN 

ADDITION, 25% MORE DENSITY OF THE EBL (D), AND A 
SYSTEMATIC SHIFT BY -18% IN THE VHE ENERGY SCALE AND, IN 
ADDITION, 25% LESS DENSITY OF THE EBL (E). NOTE THAT THE 
RESULTING SYSTEMATIC ERRORS ARE COMPARIBLE WITH THE 





STATISTICAL 


ERRORS. 








a 


^cutoff [TeV] 


A: nominal 




2.20 + 0.08 


1.44 + 0.28 


B: E+18% 




2.16 + 0.08 


1.59 + 0.29 


C: E-18% 




2.24 + 0.08 


1.26 + 0.26 


D: (E+18%) 


+ (EBL+25%) 


2.12 + 0.08 


1.61 + 0.29 


E: (E-18%) + 


(EBL-25%) 


2.20 + 0.08 


1.09 + 0.20 
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Fig. 9. — Attenuation coefficient exp(— t 77 ) for Mkn 421 
(z=0.030) using the EBL spectrum as shown in Fig. [8] 



trum and state-of-the-art cosmology (flat universe, Hub- 
ble constant H =72 km/s/Mpc, matter density f2 m =0.3, 
dark energy density Uj^=0.7) we calculated the optical 
depth r 77 for Mkn 421 . Thereby we use numeri cal in- 
tegration of Eq. 2 from lDwek fc Krennricbl ((2005) . The 



attenuation coefficients exp(—r~ n ) are shown as the func- 
tion of energy of VHE 7-rays in Fig. [5J We note that the 
attenuation coeffici ents are very similar to those from 
IStecker et all (|2006l ). 

3.3.3. The de-absorbed spectrum of Mkn 421 



9 



</> 

CM 

"E 
o 

jlio 

LU 

■o 



-9 



-10 



^10 



-11 



10 



-12 



10 



-13 



i i i i 



1" 
1 




•0 



• 


MAGIC (2004-2005) 





HEGRA (1999-2000) 


• 


HEGRA (2000-2001) 


★ 


HEGRA CT1 (2001) 


X 


CAT (1998) 




CAT (2000) 


▲ 


Whipple (2000-2001), high 


T 


Whipple (2000-2001), low 


□ 


HESS (2004) 



* I * 

r 



Q ** 



M 



j 



j 



10' 



10 ; 



J L 



10 4 

Energy [GeV] 



Fig. 10. — Differential energy spectra of Mkn 421 multip l ied by E 2 in different activity states from different experiments. The spectra 
are de-absorbed using the EBL model from Primack ct al. (20 05), upscaled by a factor of 1.5. A log-parabolic fit is performed (dashed 
lines) to determine the peak position in the SED (see Fig. 1121 . For clarity, only the highest and lowest of the published Whipple results 
are shown. Note that for the MAGIC spectrum the upper limit at 4.4 TeV is not plotted. 



The measured spectrum and the reconstructed de- 
absorbed (i.e. corrected for the effect of intergalactic 
absorption) spectrum are shown in Fig. [7j For compari- 
son reasons, the Crab Nebula spectrum is also shown. 
The de-absorbed spectrum (shown by filled black cir- 
cles) is clearly curved, its probability of being a sim- 
ple power law is 1.6 xl0~ 8 . The de-absorbed spectrum 
is fitted by a power law with an exponential cut-off: 
dN/dE = iVo(£:/0.2TeV)- a exp(-£:/i; cutoff ), a being 
the photon index, solid line in Fig.[jj The fit parameters 
are listed in the inlay of Fig.[7j The power law with a cut- 
off describes well the de- absorbed spectrum of Mkn 421, 
with a photon index a — 2.20 ± 0.08 and a cut-off en- 
ergy of E cutoS = (f -44 ± 0.28) TeV. Taking into account 
the systematic uncetrainty of 18% on the absolute energy 
scale of our measurement and in addition a guessed 25% 
uncertainty on the EBL level, we find that neither the 



photon index nor the cut-off energy substantially change 
(See Table [3]). The fitted photon index was found to be 
between 2.12 and 2.24, whereas the cut-off energy was 
found to be between 1.1 and 1.6 TeV. From this study 
we conclude that the curvature of the spectrum is source 
inherent: either at the measured flux level this cosmic 
accelerator is close to its energy limit, or there exists a 
source-intrinsic absorption. 

4. DISCUSSION 

4.1. Comparison with previous observations of 
Mkn 421 

In Fig. [TO] we show the (de-absorbed) energy den- 
sity spectrum in context with several previously pub- 
lished high statistics observations of Mkn 421. For 
a compilation of the VHE measure ments of Mkn 42 1 
we used historical data from CAT (|Piron "eTaT][200i]) . 
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Fig. 11. — Relation between the flux state at 1 TeV (deter- 
mined from a power law approximation of the spectra between 
700 GeV and 4 TeV) and the fitted photon index of published 
data as in Fig. 1101 A correlation between flux and hardness 
can be clearly seen. 
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Fig. 12. — Relation between the fitted peak position in the 
SED and the energy density at 1 TeV for published data as in 
Fig. 1101 A clear trend can be observed for the peak position to 
shift towards higher energies with increased source intensity. 



HEGRA | AliaroniaiTiO][200l [Schwcizer 20 0j|, HESS 
(lAharonian et al.| [2005). and Whipple (|Krennrich et all 
20 02) 25 . All measured spectra are de-absorbed using the 
EBL model as described in Subsection 13.3.21 The ac- 
tivity of the source during MAGIC observations clearly 
was at the lower end, and our results extend to energies 
lower than those previously observed, thus being com- 
plementary both in source intensity and energy range. 

25 For the Whipple measurements, only spectra in the high est 
and in the lowest flux state are shown in order not to clutter Fig. 1101 



All results seem consistent with each other, and all show 
significant deviations from a simple power law, which 
can not be explained by attenuation effects (the results 
are robust with respect to the EBL model within a factor 
±25%). They are, therefore, likely to be source-intrinsic. 

From the compilation of the de- absorbed Mkn 421 
spectra, it is evident that with an increasing flux state 
the spectrum becomes harder. In order to verify this, we 
fitted the spectra by a simple power law (dN/dE oc E~ Q ) 
in the overlapping energy region between 700 GeV and 
4 TeV. The resulting photon indices a as function of the 
fitted flux at 1 TeV are shown in Fig. [TTJ Evidently, 
with increasing flux the spectra hard en. Similar re- 
sults were obtained us ing Whipple data (iKrennrich et ail 
I2002T) . HE GRA data (lAharonian et all 120021 12003D. and 
CAT data (|Giebels et al.ll2006f) . 

The curvatures observed are indicative of a maxi- 
mum in energy density, and are usually interpreted 
as due to inverse- Compton (IC) scattering. The 
peak position appears to be dependent on the source 
flux intensity. We have, therefore, performed a log- 
parabolic fit for all available data. The log-parabola 
has the following parametrization: loglO(^F^) = A + 
B (log 10(E) -log 10(Ep)), with vF v = E 2 dN/dE and 
E p being the energy of the peak position. The best log- 
parabolic fits are shown in Fig. |TD] by the dashed lines. 
In Fig. [12] we compare the resulting peak positions for 
the different experiments as a function of their (fitted) 
energy density at 1 TeV. Evidently, with increasing flux 
the peak shifts to higher energy values. Future obser- 
vations at higher intensities extending to lower energies 
will have to corroborate these results. Such observations 
are part of the future MAGIC physics progamme. 

4.2. A short comment on the light curves 

In the observation period between November 2004 and 
April 2005 we observed night-to-night flux variations up 
to a factor of 2 and a maximum flux change in the entire 
set of a factor 4. No short-term flux variations well below 
1 hour, as observe d during high flaring activity in pre- 
vious experiments (jGaidos et al.l 119961 : lAharonian et al.l 
2002), were seen, although the sensitivity of MAGIC 
would allow to detect fast flares in the given flux range. 
Two equally likely explanations are that either we deal 
with large fluctuations resulting in the absence of any 
fast flare during the observation period, or fast flaring is 
a feature that occurs only when the source is very active. 
This calls for further high statistics and high sensitivity 
studies when the source is in its low flux state. 

4.3. Correlation studies 

The correlation between the 7-ray flux measured by 
MAGIC and the X-ray flux measured by RXTE/ASM 
is shown in Fig. [T31 For the MAGIC flux we take 
the nightly average above 200 GeV (see also Fig. [5]). 
For the X-ray data, we calculate the average of those 
RXTE/ASM pointings (dwells) which were taken simul- 
taneously with MAGIC data, allowing ±1 hour with re- 
spect to the MAGIC data, to increase X-ray statistics. 
Fig. [13] shows a clear correlation between X-ray and 7- 
ray data. The linear fit (solid line), forced to go through 

(0,0), has a slope of 1.4 ± 0.1 [ J ^/ £ §^], and has a X 2 
probability of 54%. The parabolic fit (dashed line) which 
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Fig. 13. — Correlation between MAGIC integral flux measure- 
ments above 200 GeV and RXTE/ASM counts for 13 nights. 
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Fig. 14. — Correlation between MAGIC integral flux above 
200 GeV and optical flux measured by the KVA telescope for 8 
nights. 



is also forced to go through (0,0) has the same \ 2 prob- 
ability of 54%. The correlation coefficient r = 0. 64^22 
(errors correspond to 1 a level) is different from zero by 
2.4 standard deviations (taking into account the non- 
linearity of errors). 

In Fig. [H the MAGIC 7-ray flux above 200 GeV is 
shown together with simultaneous KVA optical data. 
The latter have been averaged over the MAGIC inte- 
gration time. One can see a possible 7-ray/optical anti- 
correlation during the 8 nights of simultaneous observa- 
tions, however, the correlation coefficient r = — 0.591q 22 
is compatible with zero within 1.5 standard deviations. 



4.4. Comparison with models 

Given the temporal correlation between X-ray and 
7-ray fluxes, it is reasonable to infer that the VHE 
7-ray radiation is dominated by emission resulting 
from IC upscattering of the synchrotron X-ray pho- 
tons by their parent population of relativistic electrons. 
Such correlation can be modelled with a homogeneous 
synchrotron-self-Compton (SSC) model. Based on this 
model it is possible to constrain the parameter space 
of the emission region and estimate its basic parame- 
ters, the Doppler factor, D, and the rest-frame mag- 
netic field, B, of the emitting plasma in the relativis- 
tic jet. To this end we follow the p rocedure first de- 
vised bv iBednarek & Protheroel ()1997t ) for the Mkn 421 
flare of 16 May 1994, s ubsequently improved b y, e.g ., 
Tavecchio efafl (119981): IBednarek fc Protheroel fl999); 



Kataoka et all (| 19991 )7 iKatarzvnski et al.l (|2003l ). Ap- 
plication of this method requires precise simultaneous 
multiwavelength information. Since a synchrotron (X- 
ray) spectrum simultaneous with the MAGIC obser- 
vations is not available, we have to resort to previ- 
ous X-ray observations arguing that similar TeV 7-ray 
states (IC emission) should correspond to similar X- 
ray states (synchrotron emission). In fact, similar 7- 
ray spectra of Mkn 421 have already been observed 
several tim es - including the HE GRA observations in 
April 1998 {Ah aronian et al.|[l999l ) for whi ch simultane- 
ous BeppoSAX observatio ns are available (jFossati et alJ 
2000: iMassaro et al.ll200l . Here we use t he X-ray spec- 
tra an d parameterization, reported by M assaro et al.l 
(|2004h for 21 April 1998. It is also noticeable that the 
X-ray flux level between the simultaneous RXTE/ASM 
data and the BeppoSAX data used here is very similar 
(see Fig. [16]). 

The low flux state MAGIC 7-ray spectrum, reported 
here for energies at ^100 GeV, warrants a better inves- 
tigation of the crucial energy range where the IC peak 
is exp ected to occur, than in previous da ta sets. Fol- 
lowing [Bednare^&^rotheroi (|1997l . ll999i ) we then con- 
strain the allowed parameters of the emission region (D 
and B) from the ratio of the 7-ray power to the X-ray 
power, measured at their respective peak emission (see 
thick curves in the upper panels of Fig. I15|) . The ra- 
diation field density and the electron spectrum, cospa- 
tial in the blob, were derived as a function of D and 
B for a blob radius assumed equal to the light travel 
corresponding to the shortest reported variability time 
scale (for observational arguments see iTakahashi et al.l 
(2000)). We further constrain the allowed parameter 
space by arguing that the synchrotron and IC cooling 
time scales should be shorter than the observed variabil- 
ity time scale. These conditions are fulfilled above the 
dot-dashed lines (for synchrotron cooling) and on the left 
of the grey dashed line (for the IC cooling) for the 1 hr 
(upper left panel of Fig. I15P and 1 day (upper right panel 
of Fig. [T5|) variability time scales. The condition that the 
blob has to be transparent to the VHE 7-rays leads to 
a further lower bound on D by requiring that the opti- 
cal depth by pair production has to be lower than unity. 
The corresponding limits for photon energies of 100 GeV 
and 3 TeV (which define the energy range of MAGIC 
measurement) are shown in the upper panels of Fig. [15] 
as thin and thick dotted lines, respectively. One last 
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Fig. 15. — Constraints on the parameter space (Doppler factor, D, versus magnetic field strength, B) for the emission region in the jet 
of Mkn 421 based on the one-zone homogeneous SSC model. From the ratio of the 7-ray to X-ray power during the MAGIC observations 
of Mkn 421 (measured at the peaks in the synchrotron and IC spectra) the allowed values are limited to the thick curves (assumed IC 
peak at 100 GeV) or to the thin curves (assumed IC peak at 10 GeV). The left figures correspond to a 1 hour variability, whereas the right 
figures correspond to a 1 day variability. The physical conditions in the blob are limited by the electron cooling on the synchrotron and IC 
processes, optical depth for 7-rays, and maximum energy of electrons as discussed in the text. These constraints are shown by the dashed 
and dot-dashed lines. The allowed region is limited by these lines and is marked by the grey shaded area. The 7-ray spectra are calculated 
for the values of the Doppler factor and magnetic field strength marked by A, B, and C (for 1 hr variability), and D, E, and F (for 1 day 
variability). They are compared with the de-absorbed MAGIC spectrum (shown as full black circles) in the bottom figures. The fit by a 
power law with an exponential cut-off to the de-absorbed spectrum (as in Fig 01 is shown by the black solid line, whereas the measured 
spectrum is shown by the grey full squares. 



condition arises from comparing the maximum energy of 
electrons, determined by the maximum energy of syn- 
chrotron photons ~40 keV, with the maximum energy 
of the detected photons —3 TeV (see dashed line in the 
upper panels of Fig. [l3|) . These limiting conditions build 
an allowed region in the D-B plane as marked by the 
grey shaded area. The allowed parameters of the emis- 
sion region correspond to the part of the thick full curve 
inside the region limited by all these lines (see Fig. [T5|) . 
In order to determine the values of D and B more pre- 
cisely, we now calculate the 7-ray spectra for the points 
A, B, and C for 1 hr variability, and the points D, E, and 



F for 1 day variability, and compare the predicted spec- 
trum with the actual de-absorbed spectrum. From the 
lower panels of Fig. [15] it is clear that the best descrip- 
tion is provided by the blob with Doppler factor D —22 
and magnetic field B ~0.7 G (the point B) for 1 hr vari- 
ability, and D ~9 and B -0.3 G (the point D) for 1 
day variability. In order to assess how this result is sen- 
sitive on the correct energy localization of the peak in 
the 7-ray spectrum (which is in fact only limited by the 
lower energy end of the MAGIC spectrum), we show the 
allowed parameter space for the 7-ray peak at 10 GeV 
(see thin full curves in Fig. [T5|) . The constraints for the 
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Fig. 16.— The overall SED of Mkn 421 from optical wave- 
lengths through VHE 7-rays. Large symbols represent aver- 
aged data described in this paper: optical data from KVA 
(star), X-rays from RXTE/ASM (full square), de-absorbed 7- 
rays from MAGIC (full points). The grey full squa res are 
archival EGRET measurements ( Hartman et al. 19jH3). The 
grey curve in the X-rays c o rrespo nds to the log-pa rabolic 
fit ta ken fromfMassaro ct al. (2004) using BeppoSAX (Boclla 
119971 ) data of Mkn 421 taken on 21 April 1998. The two black 
curves through the 7-ray spectrum (almost indistinguishable) 
correspo nd t o the SSC model parameter sets B and D (see text 
and Fig. 1151 for details ). The grey dashed line denotes a fit 
by the SSC model as in Krawczynski et al. (2004|), see text for 
details. 



peak at 10 GeV and 100 GeV are almost the same for 
the parts of the curves inside the allowed region. It is 
interesting that the emission parameters, as estimated 
here for the low flux state of Mkn 421 (for the 1 day 
variability ti me scale), are not very differen t from those 
estimated by iBednarek fc Protherod (|1997f ) for the flar- 
ing state. This suggest that the flaring state may not 
be related to the significant change of the blob's Doppler 
factor and magnetic field strength. 

In Fig. [H we show the broadband SED of Mkn 421. 
Large symbols represent averaged data described in this 
paper: optical data from KVA (star), X-rays from 
RXTE/ASM (full square), 7-rays from MAGIC (full 
points). The grey curve in the X-rays corresponds to 
a log-parabolic fit performed by iMassaro et al.l (|2004f ) 
on BeppoSAX data of Mkn 421 taken on 21 April 1998. 
The two black curves through the 7-ray spectrum are 
almost indistinguishable and correspond to the best 
SSC model parameters for 1-hr and 1-day variability 
time scales (points B and D respectively, calcula t ed ac - 
cording to Eq. (13) in IBednarek fc Protherod (|1999D 
w ho apply the Klein-Nishina cross-section as in 2.48 
of lBlumenthal fc Gouldl (Il970l )). 

In addition, we apply the SSC code provided by 
iKrawczvnski et al"1 (|2004f ) to our dataset. The fitted 
overall SED is shown by the grey dashed line in Fig. [TH 
and the model parameters are listed in Table |4j For the 
fit, we used the simultaneous KVA, ASM and MAGIC 
data, as well as t he archival BeppoSA X observations 
from 21 April 1998 (IMassaro et alJl2004ri. In c ontra st to 
the parameters adopted in IKrawczvnski et alj (|2001h . we 
used a smaller Doppler factor (15 instead of 50), result- 



TABLE 4 

SSC MODEL PARAMETERS FOR MKN 421. THE CORRESPONDING 
PHOTON SPECTRUM IS SHOWN IN FlG. [T6l 

spherical blob with: 
Doppler factor 15 
magnetic field 0.20 Gauss 

radius of emitting region 1.6 X 10 16 cm 

injected electron spectrum: 
electron energy density 0.06 erg/cm 3 

5 < logl0(E[eV]) < 10.9 index 2.31 

10.9 < loglO(E[cV]) < 11.6 index 3.88 



ing in a somewhat larger emitting region (1.6 xlO 16 cm 
instead of 2.7xl0 15 cm), and a higher particle density 
(0.06 erg/cm 3 instead of 0.01 erg/cm 3 ). We note that 
the fitted values of magnetic field and Doppler factor 
are within the allowed range as de fined above. Remark- 
ably, the archival EGRET data (|Hartman et al.l I1999D 
matches the fit almost perfectly, suggesting an IC peak 
around 100 GeV. 

5. CONCLUDING REMARKS 

Mkn 421 was observed with the MAGIC telescope dur- 
ing several months in 2004 and 2005. Briefly, we have 
presented the following: 

• first high-sensitivity observation down to w 
100 GeV; 

• first observation of an IC peak at low flux; 

• absence of short flares below 1 hour duration de- 
spite sufficient sensitivity; 

• flux variation up to a factor 2 between consecutive 
nights and up to a factor 4 in the entire observation 
period; 

• confirmation of a source-inherent effect resulting 
in a curved spectrum after de-absorption (for rea- 
sonable assumptions concerning the EBL) in case 
of low flux intensity; 

• a strong correlation between spectral hardness 
(photon index between 700 GeV and 4TeV) and 
flux intensity, obtained by comparison of the de- 
absorbed energy spectra of various experiments 
covering different flux levels; 

• a clear trend for the peak position to shift towards 
higher energies with increased source intensity, ob- 
tained by the same comparison; 

• confirmation of a significant correlation between X- 
ray and VHE 7-ray intensity during a state of low 
to medium intensity; 

• a hint that different flaring states result from differ- 
ences in electron populations (electron spectrum) 
rather than from significant change of the blob's 
Doppler factor and magnetic field strength. 

We add the following conclusions. The flux state 
was found to be comparatively low, ranging in inten- 
sity between 0.5 and 2 Crab units when integrated above 
200 GeV. While clear night to night variations were 
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found, the intra-night light curve, binned in 10-minute 
time intervals, does not show significant variations, al- 
though several nights are only marginally compatible 
with a constant flux. They do not show a discernible 
structure, though, and seem not associated to an overall 
flux different from that of perfectly quiescent nights. We 
note that MAGIC is sensitive enough to detect variabili- 
ties on the 10-minute time scale at such a low flux level. 
A clear correlation (r = 0.64j^ 2 |) between X-rays and 
7-rays was found, while no correlation was seen between 
optical and 7-rays. This supports a leptonic origin of 
the 7-rays from Mkn 421. The energy spectrum result- 
ing from the combined MAGIC data, corrected for the 
extragalactic absorption, suggests the presence of an IC 
peak at about 100 GeV. The spectrum is clearly curved 
at energies above 1 TeV, and can be fitted by a power-law 
with an exponential cut-off. The overall SED observed 
in the observed flux state can be well described by a 
homogeneous SSC model provided that the emission re- 



gion moves with a Doppler factor ^9 and its magnetic 
field strength is ~0.3 G for a 1-day variability time scale. 
Surprisingly, these parameters do not differ substantially 
from those estimated for the emission region of Mk n 421 
during a strong flare (jBednarek &: Pro thcroc 1997) . The 
fit wit h an alternative SSC code of iKrawczvnski et al.l 
(|2001h lead to similar Doppler factor and magnetic field 
values. 
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